Star copolymers in porous environments: scaling and its manifestations 
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We consider star polymers, consisting of two different polymer species, in a solvent subject to 
quenched correlated structural obstacles. We assume that the disorder is correlated with a power- 
law decay of the pair correlation function g(x) ~ x~ a . Applying the field-theoretical renormalization 
group approach in d dimensions, we analyze different scenarios of scaling behavior working to first 
order of a double e — 4 — d, <5 = 4 — a expansion. We discuss the influence of the correlated disorder 
CN) ' on the resulting scaling laws and possible manifestations such as diffusion controlled reactions in 

the vicinity of absorbing traps placed on polymers as well as the effective short-distance interaction 
O ■ between star copolymers. 
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PACS numbers: 82.35.Jk, 36.20.Fz, 64.60.ae, 64.60.F- 



I. INTRODUCTION 



Understanding the behavior of polymer macromolecules in solutions in the presence of structural obstacles is of 
great interest in polymer physics. The presence of defects often leads to a large spatial inhomogeneity and may create 
pore spaces of fractal structure Such situations can be encountered in studying, e.g., polymer diffusion through 
I microporous membranes Q or within colloidal solutions Q. 

Solutions of polymer macromolecules in disordered environment are subject to intensive studies. Numerous sim- 
ulations jJQ and analytical studies 0, HoHT^ I have focussed on the case of uncorrelated structural defects at the 
percolation threshold of the remaining accessible sites, this situation is shown to alter significantly the universal be- 
havior of polymer macromolecules. Recently, another special type of disorder has been brought to attention, which 
display correlations in mesoscopic scale. This case can be described within the frames of a model with long-range- 
correlated quenched defects, considered in Refs. fl7l - fl9j in the context of magnetic phase transitions. Here, structural 
defects are characterized by a pair correlation function g{x), which in d dimensions falls off at large distance x 
C^Q ■ according to a power law: 

t-H ■ 

<N ; g(x) ~ x~\ (1) 

O 

In general, any value of < a < d can be realized by defects that form clusters of fractal dimension df = d — a. 
For integer dimension df these include the following special cases: uncorrelated point-like defects (df = 0), mutually 
uncorrelated straight lines of random orientation (df = 1), mutually uncorrelated planes of random orientation 
(df = 2). The influence of such long-range correlated defects on the universal properties of single polymer has been 
analyzed within the renormalization group approach in Refs. (20L l2l|. 

To describe the universal properties of polymer chains in good solvents, one may due to universality in the long 
chain limit consider the model of self-avoiding walks (SAWs) on a regular lattice [22], Hit- I n particular, the average 
square end-to-end distance (R 2 ) and the number of configurations Zn of SAWs with N steps obey in the asymptotic 
limit N — » 00 the following scaling laws: 

(Rl) ~ N 2v , Z N ~ R 2 -^ 1 ^. (2) 



Here, the second equation shows the power law in terms of the effective polymer size R = */ (R 2 ) ~ N", v and r\ are 
universal exponents that only depend on the space dimensionality d. For d = 3, high order renormalization group 
estimates are [H v = 0.5882 ± 0.0011 and 77 = 0.0284 ± 0.0025. 

The theory can be generalized to describe star polymers, which consist of / linear polymer chains or SAWs, linked 
together at their end-points. The study of star polymers is of great interest since they serve as building blocks of 
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polymer networks [2|| [26[ and can be used to model complex polymer micellar systems and gels (27l - l29| . For a single 
star with / arms of N steps (monomers) each, the number of possible configurations obeys a power law in terms of 
the size R of the isolated chain of N monomers [2a, : 



R ri f -fm, (3) 



Here, the exponents rjf are universal star exponents, depending on the number of arms / (771 = 0, 779 = \ jv — 2 — rf). 
Scaling properties of star polymers are well studied both numerically [3014351 ] and analytically l36l - l42T | . It has been 
shown that the presence of long-range-correlated disorder may have interesting consequences for the scaling properties 
of polymer stars such as entropic separation of polymers according to their architecture [431 ] . 

Linkin g togeth er polymers of different species, we receive non-homogeneous star polymers with a much richer scaling 
behavior 44j-47j. A particular case is the star copolymer, consisting of polymer chains of two different species. It has 
been shown [45], that the number of configurations Z^f 2 of a copolymer star with fx arms of species 1 and arms 
of species 2 scales as: 

Z flf2 ~ (i?)'7/l/ 2 -/m20-/27702 ) ( 4 ) 

where rjf 1 f 2 constitutes a family of copolymer star exponents. These exponents are universal and depend only on 
space dimension d and the number of chains /1, /2, as well as three different types of fixed points that govern the rich 
scaling behavior [44| . 

Depending on the temperature, a situation may occur, where one or more of the inter- or intrachain interactions 
vanishes. Indeed, for each polymer system one finds a so-called 0-temperature, at which attractive and repulsive 
interactions between monomers compensate each other (see, e.g., [22|, Hj]). Such a polymer chain can effectively be 
described by a simple random walk (RW). In this case, scaling laws (Eq. ©) hold with exponents: v = 1/2, r) = 0. 
As a result, for example, there may be only mutual excluded-volume interactions between chains of different species, 
while chains of the same species can freely intersect. That is, some species behave effectively like RWs. Within 
a copolymer star, the subset of chains of such species builds up a substar of random walks, possibly avoiding the 
second part of the star, which can be either of random walks or self-avoiding random walks (see Fig. [T]). Cates and 
Witten [48{ have shown, that this situation can also be interpreted as describing the absorption of diffusive particles 
on polymers. 

Another example, where star exponents gove rn physical behavior concerns the short-range interaction between 
cores of star polymers in a good solvent [251 l49l |50| . The mean force Ff g (r) acting on the centers of two stars with 
functionalities / and g is inversely proportional to the distance r between their cores: 

Ffg (r)=k B T^L, (5) 

with ksT denoting the thermal energy, Qf g the universal contact exponent, related to the family of exponents of star 
polymers by scaling relations: 

®fg = Vf + Vg-rif+g- (6) 



FIG. 1: (color online) Schematic representation of copolymer stars consisting of two polymer species (denoted as red and blue). 
Solid lines present species behaving like SAWs, dashed lines present RWs. The two different sets in each example may further 
be either mutually avoiding or mutually "transparent". 
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We are interested to generalize this relation to the case of copolymers, and to analyze the impact of disorder on 
mutual interactions between two star copolymers. 

The questions of the influence of correlated disorder in the environment on the scaling behavior of star copolymers 
and the resulting consequences remain so far unresolved and are the subject of the present study. We will also analyze 
the spectrum of scaling exponents in particular for the above mentioned process of trapping diffusive particles in the 
vicinity of absorbing polymers in disordered environments. 

The paper is organized as follows: in the next section we will give a field-theoretical representation of the model 
Lagrangean. The field-theoretical renormalization group method, which we use to find the qualitative characteristics 
of scaling behavior, is shortly described in Section 3. In Section 4 we discuss the results obtained. We finish by giving 
conclusions and an outlook. 



II. THE MODEL 



Let us consider a polymer star with / arms of different species in a solvent. We are working within the Edwards 
continuous chain model [HLH^, representing each chain by a path ri(s), parameterized by < s < Si, i = 1, 2, . . . , /. 
The central branching point of the star is fixed at ri(0). The partition function of the system is then defined by the 
path integral [26| : 

Zf(Si) = I D[n,...,r f ] 

f 

xexp n^MO) -^(0)). (7) 

Here, a multiple path integral is performed for the paths ri,...,r/, the product of (^-functions reflects the star- 
like configuration of / chains, each starting at the point ^i(O), T-Lf is the Hamiltonian, describing the system of / 
disconnected polymer chains: 

i=l 



1 s r 

- £ u%j dr Pi (r) Pj (r), (8) 



i<i=i 

r-S.. 



where Pi{r) = J Q ' ds 5 d (r — rj(s)) and u®j is a symmetric matrix of bare excluded- volume interactions between chains 
i and j. 

The continuous chain model ([7]) can be mapped onto a corresponding field theory by a Laplace transform in the 
Gaussian surface Si to the conjugated chemical potential variable (mass) jji [ID, |44j : 

%(£0 = / II d ^ expH^-Z/^)- (9) 

J b 

One may then show that the Hamiltonian H is related to an m-component field theory with a Lagrangean C in the 
limit in — y 0: 

/ 



i f r 

+ 1 L <j d^ 2 (*V 2 (x), (10) 



»<i=i 



where ipf 1 = {(fj, . . ■ , ip™} and \n are bare critical masses. On the base of the Lagrangean (Tl7J|) the one-particle 
irreducible vertex functions r^ L -* of the theory can be obtained: 

8(£) r Z--M = J e iqiri dr 1 ---dr L (<p il (r 1 )...<p iL (r L ))f PI , (II) 
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where only those contributions, that have non-vanishing tensor factors in the limit m —y are kept. 

The Laplace-transformed partition function Zf(p,i) has a vertex part, which is defined by the insertion of the 
composite operator Yii fi- 

S(p + J2<lj)^ f (p,<lU---,Qf)= (12) 



Apro+m) drQdri dr/ ^ l(ro) . . . <p f ( ro )<pi(ri) ■ ■ ■ <p f (r f )) 



c 



Let us note that we are interested in the case of a copolymer star, having j\ chains of one species and fi of another, 
so that f\ + $i = f. To keep notations simple we will consider in the following discussion only two fields tpi and if2, 
corresponding to two different "species". Thus, in (|10l) we have interactions Uii,U22 between the fields of the same 
"species" and u±2 between different fields. The composite operator in (fT2")l has the form of a product ((fii)^ 1 (ip2)^ 2 ■ 

We introduce disorder into the model (jlOl) . by redefining /if — > /if + Sfli(x), where the local fluctuations Sfii(x) 
obey: 

((6fjLi(x)))=0, 

((5ni(x)5(j. j {y)}) = gij(\x - y\). 

Here, ((• • •)) denotes the average over spatially homogeneous and isotropic quenched disorder. The form of the pair 
correlation function g(x) is chosen to decay with distance according to the power law ([l}. 

In order to average the free energy over different configurations of the quenched disorder we apply the replica 
method to construct an effective Lagrangean: 

2 



A 1 1 



i=l a— 1 
2 n 

+ E Ei^f(^) 2 as) 

i<j — l a— 1 
„ 2 n 

dxdy E 9iA\x-y\m?i$?- 

J i<i=l a,/8=l 

Here, the coupling of the replicas is given by the correlation function g(x) of Eq. ([I]), Greek indices denote replicas 
and the replica limit n —y is implied. 

For small k, the Fourier-transform gij(k) of gij(x) ([T]) reads: 

^•(fe)~^+<.|fcr d . (w) 

Thus, rewriting Eq. (|13[) in momentum space, one obtains an effective Lagrangean with 9 bare couplings: w^, uS^j 
u 12j ^22) u i2) w ii > ^22; u 'i2- As ^ was pointed out in Ref. [T3| ■ once the limit m, n — >■ has been taken, the u®j 
and terms acquire the same symmetry, and an effective Lagrangean with couplings (u°- — = u° ) of 0(mn = 0) 
symmetry results. This leads to the conclusion that weak quenched uncorrelated disorder i.e. the case a = d is 
irrelevant for polymers. Taking this into account, we end up with only 6 couplings in an effective Lagrangean: u^, 
u 22i u i2> ^li) w 22> w i2- F° r a < d, the momentum-dependent coupling w { -jk a ~ d has to be taken into account. Note 
that gij(k) must be positively definite being the Fourier image of the correlation function. Thus, we have w?- > 
for small k. Note, that the couplings it^ should be positive, otherwise the pure system would undergo a 1st order 
transition. 



The resulting Lagrangean in momentum space then reads: 



n 2 

^« = ^EEE[ fc2 +^] «(fc)) 2 



a — 1 2—1 k 
2 / n 

■ E E if E 6 (ki+k2+k 3 +k 4 ) tf(k 2 ) 0f(k 3 ) ^(fc 4 ) (15) 

i<j=\k\h^ \ a — 1 



£ |fcl+fc 2 | a -^(*:i+* a -fc3-fc4)^*(*l)^(*2)^(*3)^(*4 
a,P=l 
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In the next section, we apply the held-theoretical renormalization group approach in order to extract the scaling 
behavior of the model (fT5|) . 

III. RENORMALIZATION GROUP APPROACH 

We apply the renormalization group (RG) method [53[ in the massive scheme renormalizing the one-particle ir- 
reducible vertex functions, in particular r^ 2 \r^ and T 2 ' 1 , as well as the vertex function T*^ 1 '^, with a single 
(pi)-' 1 (<P2) insertion. Note that the polymer limit of a zero component field leads to an essential simplification: each 
field (fi, mass /ii and coupling renormalizes as if the other fields were absent. The renormalized couplings Uij,Wij 
are given by: 

u% = t j, e Z^Z ii un ) i = l,2 (16) 
wl = ii s Z-?Z u w iU i = l,2 (17) 
u° 12 = fi £ Z-^Z-^Z 12 u 12 , (18) 
<2 = V S Z-'Z-iZ 12 w 12 . (19) 

Here, [i is a scale parameter, equal to the renormalized mass, and parameters £ = 4 — d, <5 = 4 — a. The renormalization 
factors Z have the form of power series, the coefficients of which are calculated perturbatively order by order. 
The star vertex function r*" 1 '-^ 2 ^ is renormalized by a factor Z„f 1 j 2 : 

zMzWz.mTM*) = M (A+/ 2 )(e/2-D+4- e . (20) 

The variation of the coupling constants under renormalization defines a flow in parametric space, governed by 
corresponding /3-functions: 

Pun (uij,Wij) = fJ>—Uij, [3 Wi] (uij,Wij) = fi-^-Wij, i,j = 1, 2. (21) 

The fixed points (FPs) of the RG transformation are given by the solution of the system of equations: 

«j, <j) = 0, fi Wij (u*j, w^) = 0, i,j = 1, 2. (22) 

The stable FP, corresponding to the critical point of the system, is defined as the fixed point where the stability 
matrix possesses eigenvalues {Aj} with positive real parts. 

The flow of the renormalizing factors Z in turn gives rise to RG functions r\ lpi and ?7*/ 1 / 2 as follows: 

li-^-lnZ Vz = Tj^UijjWij), (23) 

/^ ln ^*/i/2 = V*fih(uij,w i:j ). (24) 

At the FP of the renormalization group transformation, the function r] ipi describes the pair correlation function critical 
exponent, while the functions Ty*/^ define the set of exponents for copolymer stars: 

V = V<Pi(uZj,Wy) (25) 

Vfih=V*hh{u-ij,Wij)- (26) 

In the next section, we will present expressions for the /3 and rj functions, together with a study of the RG flow and 
the fixed points of the theory. 

IV. THE RESULTS 

A. Fixed points and scaling exponents 

According to the renormalisation group prescriptions, we obtain the RG functions of the model (|15[) within a 
massive scheme up to the one-loop approximation: 

\(26-e)wll 3 , (27) 





4 2r" 






Pun = - £ 




- 52u li w il 


h + \h 
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+ e-wuUiih, i = 1,2; 



(28) 



2 1 

"12 - ^u\ 2 Il - -Ul 2 {un + U22Kl 



u 12 w X2 l 2 + 5^12 (wii + w 2 2)h + ■xuuiwn + w 22 )I 4 ) 



+(25 - e) 

1 



1 2 r 1 / N 

-w 12 / 3 + -W 12 {Wii + w 2 ) 
3 6 



1 2 r 1 2 r 

">12 + ~U>l 2 -<2 + ^W u h 



+e 



W l2 Ul2h + \wi2{uu + U 22 )Ii + ^Wl2(^ll + Wi 2 )I 2 ) 

3 b b 



(29) 



(30) 



Note, that expressions for /3 Uii ,f3 Wii restore the corresponding RG functions for a single polymer chain in long-range 
correlated disorder 0, [2l|. Here, Ii are the loop- integrals: 

dq 



h 
h 
h 



d 



(g 2 + l) 2 ' 

dqq a ~ d 
(g 2 + l) 2 ' 

{q 2 + l) 2 ' 

dqq a - d 



[q + k] 2 + l) 



(31) 



k 2 =0 



We make the couplings dimensionless by redefining Uij — Uijfi and Wij — Wij/i a ~ , therefore, the loop integrals 
do not explicitly contain the mass. Besides, we absorb geometrical factors Sd, resulting from angular integration into 
the couplings. 

Additionally, we need the RG function Tj„ l f 1 f 2 (uij,Wij), which we find in the form: 

V*hf3 =-e( u n g h+u 22 - h +u u — — h 

\ 6 J 2+W22 g h + wi 2 — — 1 2 J . (32) 

The perturbative expansions for RG functions may be analyzed by two complementary approaches: either by 
exploiting a double expansion in the parameters e = 4 — d, 8 = A — a [l?], H3, EH , or by fixing the values of the 
parameters d, a [20 ] . Let us note, that within the one- loop approximation the latter method can not give reliable 
results [2(3], and we exploit the double expansion in e = 4 — d, 5 — 4 — a for a qualitative analysis. The resulting 
expressions for /3- and 77-functions read: 



-sua + ^u 2 - 2u u w u + ^u> 2 



SWi; 



2 2 

2 W ti + 3 M « W ": i = l,2; 



2 2 1 , x 1 , 

-£Mi2 + -M 12 + -"l2("ll + U22) - 2^12(^11 + w 22> 



1 2 1 

-M12W12 + -w 12 + -Wl2(Wn + W22), 
3 

1 2 1 

Pw 12 = -OW12 - ~W 12 + ~U 12 W 12 

+ \w X2 (un + u 22 ) - \wi 2 {wn + w 22 ), 
6 6 

/i(A-i) , x Mh-i) ( , 

V*flf 2 = g ("11 - Wll) g ("22 ~ W 22 ) 

, /2(/2-l) 
^ "12 H ^ Wl2- 



(33) 
(34) 

(35) 

(36) 

(37) 
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Substituting Eqs. (I33|) - (l36l) into (|22|) . we find a number of fixed points, corresponding to different scenarios of the 
scaling behavior of the model. 

Pure solution First, let us consider the case, when disorder is absent (ton = im? — w i2 — 0) and we recover the 
problem of the so-called ternary solution of two polymer species in agood solvent [44j . Solving the equations f3 Uij = 0, 
i,j = 1, 2, we find eight fixed points in correspondence with Refs. [45144 7| . The trivial FPs: Gofu*! = u'22 = u\ 2 = 0), 
t/o(u*j ^ 0, 1*22 = "12 = 0), U (u22 7^ 0, u\ x — u* 2 = 0) and So( u ii = u 22 7^ 0' u \2 = 0) describe sets of two mutually 
non-interacting polymer species. More interesting are the FPs denoted as G, U, U , S, describing two mutually 
interacting species, their coordinates are given in the upper part of Table HI Corresponding values of the exponents 
Vhh read: 

a -(hh)e 
Vhh ~ 2 ' 

u u' -A(/i+3/ 2 -l)e 

^hh ~ Vf 2 ,h ~ § ' 

S -(/l+/2)(/l+/2-l)g ,„ R , 

Vhh ~ g • ^1 

Note, that Vfafe just recovers the exponent of a homogeneous polymer star with / = fx + f% arms. The values of 
these exponents are known up to 4th order of the e-expansion (2(| H3] and in the fixed d approach [45[ . 

Solution in the presence of long-range correlated disorder. Next, let us turn on the disorder. Apart from the eight 
FPs listed above, now we have a whole set of new FPs describing two polymer species in the case, when one or both 
of the species feel the presence of long-range correlated disorder. Indeed, to find these FPs one has to solve the 
system of 6 second-order equations (|22"T) with the /?- functions given by (|33]l — (|36| . In principle, this may lead to 2 6 
solutions [55l ] . In what follows we consider only four nontrivial points, corresponding to copolymer stars of mutually 
interacting species, both feeling the presence of disorder, which are of foremost interest (see Table |T]). These FPs 
describe particular situations of two mutually interacting sets of RWs (Gl), SAWs (Sl) and two interacting sets of 
RWs and SAWs (Ul, U l ). Note that due to the special form of the /3-functions the fixed points with u* { = 0, w* { ^ 
do not exist, i.e. one cannot describe simple random walks in the media with long-range-correlated disorder. 

We are interested in the points, which are stable in all coordinate directions. After analyzing the stability and 
physical accessibility of all the points, we come to the conclusion, that only the FPs 5* and Sl are stable in all 
directions and their stabilities are determined by the conditions: 

• fixed point S is stable for e > 25, 

• fixed point Sl is stable for 6 < e < 28. 

Although the remaining FPs (Gl, Ul and U L from the Table 1) are unstable, they can be reached for 5 < e < 26 
under specific initial conditions. In particular, Gl is reachable from the initial condition u\x = U22 — wn = W22 = 0, 
Ul is reachable for u 2 2 = W22 = and U L for Un = u>n = 0. Substitution of these FPs values into the expansion 
([57]) results in the following estimates for r]f 1 f 2 : 

Vh% = -(AW 

u L _ u' L _ -A(/i+3/ 2 -l)£ 

Vfih - 7 '/2,/i _ A 



TABLE I: Non-trivial fixed points of the model (fT51) . 

Ml! U22 ^12 ^11 U'22 W12 

G ^ 

u f f 

u' % f 

s f f f 

Gl <gy ^ffi 

rr 3S 2 n 9<S 2 3S(e-2S) n 9<5(2<5-e) 

UL 2(e-S) U 4(e-<5) 2(<5-e) U 4(e-S) 

rr' n ' Ag2 9s2 n 3^(e-2^) os{2s-e) 

U L u 2( £ -«) 4(e-i) U 2(«-e) 4(e-«) 

c 3S 2 3S 2 3S 2 36(e-2S) 36(?-26) 3S( s -26) 

° L 2{e~5) 2(e-5) 2(e-5) 2{S-e) 2(S-e) 2{S-e) 
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Here, f\£f % gives the exponent for the homogeneous star with /1 + /2 arms in solution in long-range-correlated disorder, 

^l\h anc ^ ^Jih Ascribe f 2 random walks, interacting with f\ RWs and with /1 SAWs respectively, in long-range- 
correlatcd disorder. All this leads to a variety of new scaling behavior for copolymer stars in a disordered medium. 



B. Diffusion-limited reaction rates 



Let us consider the /i-arm star polymer with arms of linear size R s and absorbing sites all along these arms. At the 
center of the star a particular absorbing trap is placed. Free particles A which diffuse in solution are trapped or react 
at these sites. We are interested in the reaction rate fc/ 2 of simultaneously trapping fi randomly walking particles A. 
This rate is proportional to the averaged moments of the concentration p of the particles near this trap and scales as 

k h ~ {pt>) ~ R s Xm \ (40) 

This process is an example of a so-called diffusion-limited reaction H3 > with the rate depending on the sum of the 
diffusion coefficients of the reactants [58| . As far as the presence of disorder lowers the diffusion coefficients [HI, [6(| , 
it is predicted to lower rates of association in diffusion-limited circumstances. It is interesting to check this prediction 
analytically, analyzing the behavior of star copolymers in long-range correlated disorder. In terms of the path integral 
solution of the diffusion equation, one finds that to calculate the rate of a reaction at the absorber that involves f 2 
particles simultaneously one needs to consider f 2 RWs that end at this point. The moments of concentration in Eq. 
([4TIf are thus defined by a partition function of a star comprising f 2 RWs [2^, [26[ . Finally, introducing the mutual 
avoidance conditions between the absorbing star and a "star" of diffusive particles one ends up with the problem of 
calculating the partition function of a copolymer star with two species /1, f 2 . By means of the short-chain expansion 
[49l | the set of exponents r]f 1 f 2 in © can be related to the exponents \f±f 2 in (jlO)) [45U471 . l6l| : 

A /l/2 ''/l/2' 

Based on these relations, the resulting values for the pure solution read [6l| : 



y RW pur e _ s f, 

A /i/ 2 ^ 2 /l/2 ' 



aSSST- = jhh. (42) 

Let us note, that the case fx = 2 corresponds to a trap located on the chain polymer, whereas /1 = 1 corresponds to 
a trap attached at the polymer extremity. 

Corresponding values for the exponents defining these processes in an environment with long-range correlated 
disorder can be obtained by substituting Eqs. (f39| into (|4T|) : 



Comparing relations P^|) and P5|) at fixed values e = 1 (d = 3) and varying the parameter 6, one notes that the 
presence of correlated disorder results in an increase of the exponents A. Moreover, the stronger the correlation of 
defects, the larger is A. Recalling the definition (T4TH) . we immediately conclude that, as expected, the presence of 
long-range correlated disorder results in lowering the rates of diffusion-limited reactions. The crucial point is that 
while long-range correlated disorder apparently does not influence the RW itself (there is no new fixed point with 
uu = 0, Wu =/= 0), the fact that the absorbing polymer changes its conformation and fractal dimension in the LR 
background leads to a change of the diffusive behavior of particles being absorbed (or catalyzed) on the polymer. 

Let us analyze several particular cases: 

• For a given /i-star absorber i.e. a reactive site placed at one end of an otherwise absorbing polymer increasing 
the size R s by a factor of / changes the reaction rate to k'^^ ~ (Zi? s ) — A/ i / 2 , so that: 

k 'hfJ k hh~ rX ' 112 - ( 44 ) 
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Increasing the size of the polymer thus leads to a reaction rate decrease by a factor of I *fita . Since Xf ± f a 
is larger than Xf"j 2 , we conclude, that the presence of long-range correlated defects makes the reaction rate 
decreases more slowly as compared to the pure solution case. 

• For a fixed number f% of particles to be trapped simultaneously the effect of attaching f 1 additional arms to an 
/i-arm star absorber decreases the reaction rate: 

k f 1+ fUJ k ^~ R ^ X,1+/lf2 > ^ 

as far as f 2 > ^/i/ 2 - This decrease is suppressed to some extent in the presence of long-range correlated 

defects. 

• For a given /]-star absorber an increase of the number of particles to be trapped simultaneously results in a 
decrease of the reaction rate: 

since A,, + v > ^f 1 f 2 - Again, presence of disorder makes the reaction rate decrease more slowly as compared 
to pure case. 



C. Interaction between star copolymers 



The effective interaction between two star copolymers at short distance r between their cores can be estimated 
following the scheme of Refs. [H, Hi|[5(|, based on short distance expansion. The partition sum Zf 1 f 2 9lff2 (r) of the 
two stars with f = fx + f% and g — g\ + gi arms of species 1 and 2 at small center-to-center distances r factorizes into 
a function C/ x / 2 gi92 (r) and the partition function Zf 1+gi f 2 + g2 of the star with /i + g\ arms of species 1 and + g2 
arms of species 2 which is formed when the cores of the two stars coincide: 

Z fl h 91 92 ( r ) - C/l h 91 92 ( r ) Z h +91 h +92 ' ( 47 ) 

For the function Cf 1 f 2 gig2 {r) it was shown [25L l49j that power-law scaling for small r holds in the form: 

G/l/2 9l92M-r 9/l/ ™- (48) 

To find the scaling relation for this power law, we take into account Q and change the length scale in an invariant 
way by: r — > £r, R —> £R. Eq. (|47l) then can be written as: 



£ fif 2 gis 2 £ r lfi+ai /a+92 _ (/i+S r i)'»D-(/a+ffa)»7oa^ f 2 +g 2 - (49) 



Collecting powers of i provides the scaling relation for the contact exponent: 



e /i/ 2 9192 = I/1/2 - hmo - M02 + Vgig 2 ~ 9imo - 32??02 - (50) 
(Vh+91 / 2 +92 - (/1 + 51)^20 - (/a + 92)^02) = 
Vhf2 ' r \g\g2 — Vfi+gi h+g2- 
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For two star copolymers at a distance r between their centers the mean force Ff 1 f 2 gig2 (r) acting on the centers can 
be derived as the gradient of the effective potential V e ff(r) — ~k J T\og[Zf 1 f 2 giS2 (r) / \Z f 191 Z / 292 )]. For the force at 
short distances r this results in (62[: 

F flhgi92 (r) = kjT e ^ g \ (51) 

For two mutually interacting star copolymers we have the three following nontrivial situations. First one may have 
two stars, each consisting of two species (with numbers of arms /i, f 2 and gi, g 2 respectively) all behaving as mutually 
avoiding SAWs (see Fig. [5^). This situation is equivalent to two SAW star polymers of the same species. A second 
possible situation is the interaction between two star copolymers, the first containing f\ SAWs and f 2 RWs, another 
<7i SAWs and g 2 RWs (Fig. [2p). Thirdly, one may have two stars, each consisting of two species (with fi, / 2 and gi, 
g 2 arms respectively) all behaving like RWs but with mutual avoidance between the species (Fig. [5t). It is easy to 
check, that any other case can be represented in terms of these three nontrivial examples. E.g., putting f 2 — in the 
case corresponding to Fig. ^jp, we obtain a homogeneous /i-arm star polymer interacting with a star copolymer, etc. 

Taking into account Eqs. (j38| . (|39|) we find the following contact exponents corresponding to the three nontrivial 
situations described above. 

1) Pure Solution 

7 ffi+f 2 + V gi+g2 ~ Vf 1+ f 2+gi+g2 = + /a)(ffi + 92) (52) 

Vf lh + < 92 - €+91 h+ 9 2 = f (2/i5i + 3/i32 + 3 fl i/ a ) (53) 
7 ?/i/ 2 + vf l9a ~ Vh+ai /2+92 = f (A02 + h9x)- (54) 

4 L +h + - = 2 (/l + h){gi + 92) (55) 

<fhh + € L 92 ~ rt +gi h+92 = |(2/i5i + 3/xft, + 3. 9l / 2 ) (56) 
€h + <%, h+92 = S ^2 + hax). (57) 

Qualitative estimates for the contact exponents in d = 3 can be found by direct substitution of e — 1 in the above 
relations. To discuss the physical interpretation of these results, let us consider Fig. |3l comparing the cases of pure 
lattice and LR disorder with a — 2.2 and a = 2.7. Fig. [3^ presents the contact exponent £ 0g2 governing the 
interaction between a star copolymer and a homogeneous star with g 2 arms of RWs. We fix g 2 — 8 and change /1 
and f 2 in such a way that f\ + f 2 = 8. The case f\ = 0, describing two stars of RWs, results in zero value contact 
exponents and thus the absence of interaction. Increasing the parameter fx (SAW component) leads to the gradual 
increase of the strength of the interaction. For j\ = 8, we have a star of SAWs interacting with a star of RWs with 
maximal interaction strength. Fig. depicts the situation of a star copolymer interacting with a star of g\ = 8 
SAWs. Again, we change j\ and f 2 as above. The case fi — describes a star of SAW interacting with a star of RWs 
and is a particular case of Fig. [3^ described above. Increasing /1 leads to a gradual decrease of the strength of the 
interaction. For fi = 8, we have two interacting stars of SAWs with minimal interaction strength. Fig. [3J; depicts a 
situation of two interacting star copolymers with /1 , f 2 and g\ , g 2 arms correspondingly. We fixed g\ = g 2 = 4 and 
again change /1 and f 2 as described above. The case fi — describes a copolymer star interacting with a star of 
RWs. Increasing the parameter /1 leads to a gradual increase of the strength of the interaction, until it reaches its 
maximal value at fi = 8, corresponding to the interaction between a star copolymer and a star of SAWs. The case 
fi = f 2 = 4 describes the interaction between two identical copolymer stars. 

Finally, we conclude, that in all situations considered above, the presence of correlated disorder leads to an increase 
of the contact exponent. The stronger the correlation (the smaller the value of correlation parameter a), the stronger 

( S S'} 

is the interaction between polymers in such an environment. Let us recall, that the exponent Sl92 corresponds 

to the situations of two interacting homogeneous polymer stars of / = /1 + f 2 and g = g\ + g 2 arms in solution in the 
presence of long-range-correlated disorder. This problem has previously been analyzed [43] using a two-loop expansion 

( S 

series for 9j 9 in d ~ 3. The quantitative estimates obtained predict a decrease of the contact exponents with the 
strength of the disorder correlations in contrast to our present e, ^-expansion results. Revising the resummation as 
performed in [43| we conclude that the number of terms in the two-loop expansion is probably too small to rely on 
those quantitative results. 



/1/2 9192 
/1/2 91 92 
hh 9l 92 



2) Presence of LR disorder 



Q(S S)l 
/1/2 Si 92 

/1/2 5152 

/1/2 3132 



11 



FIG. 3: Contact exponents of interaction between a copolymer star with /i SAWs and fa RWs and: an 8-armed star of RWs(a); 
an 8-armed star of SAWs(b); a copolymer star with 4 arms of SAWs and 4 arms of RWs in d = 3. Boxes: pure case (a = 3), 
circles: a — 2.7, triangles: a = 2.2 
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CONCLUSIONS 



In the present paper, we have studied the scaling properties of copolymer stars, consisting of f± arms of polymer 
species 1 and fa arms of species 2 in a solution in which one or more of the intra- and interspecies interactions are 
found to be at their 0-point with the further complication of a disordered environment with correlated structural 
defects. We assume that the disorder is correlated with a power-law decay of the pair correlation function g{x) ~ x~ a 
at large distance x. This type of disorder is known to be relevant for simple polymer chains [20I. f2lj| and homogeneous 
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polymer stars [43[, and we address the question of the scaling of copolymer stars in this situation. 

Considering the /i-arm absorbing star polymer with a special trap placed at the center of the star where /2 free 
particles (RWs) are to be trapped simultaneously, the reaction rate of this diffusion-limited reaction is found to scale 
with exponents, connected to the spectrum of critical exponents rjf 1 f 2 of star copolymers [48j. Such a process is 
an example of a so-called diffusion-limited reaction, with the rate depending on the sum of diffusion coefficients of 
the reactants. Another example, where star exponents govern physical behavior concerns the short-range interaction 
between the cores of star polymers in a good solvent. The present study aims to analyze the impact of structural 
disorder in the environment on these processes. 

In the frames of the field-theoretical renormalization group approach, we obtain estimates for the critical exponents 
rjf 1 f 2 up to the first order of an e = 4 — d, 6 = 4 — a-cxpansion, which belong to a new universality class. In particular, 
this enables us to conclude that the rates of diffusion-limited reactions are slowed down by the presence of long-range- 
correlated disorder. The crucial point is that while long-range correlated disorder apparently does not influence the 
RWs and thus the universal behavior of diffusion itself, the fact that the absorbing polymer changes its conformation 
in the LR background leads to a change of the rate with which particles are absorbed (or catalyzed) on specific sites 
of the polymer. 

The contact exponents, governing the repulsive interaction between two star copolymers in correlated disorder, are 
found to be larger than in the pure solution case. The stronger the correlation of the defects, the stronger is the 
interaction between polymers in such a disordered environment. 
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